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ABSTRACT 

In this paper we investigate the environment and reionization process around the highest redshift 
QSOs having Gunn-Peterson troughs (z > 6.1). Starting with the cosmic density perturbation and 
structure formation theory and the fact that the highest redshift QSOs are located in rare overdense 
regions, we show that the halo formation, gas distribution, and star formation around QSOs are biased 
from those of the cosmic average. We argue that a significant fraction of hydrogen in the Stromgren 
sphere around QSOs is ionized by photons from stars and that only about several percent to at most 
10%— 20% of the total hydrogen is left (e.g., in minihalos, halos, or high-density subregions) to be 
ionized by QSO photons. The cosmic average neutral hydrogen fraction at z ~ 6.2 — 6.4 should also 
be smaller than the upper limit of 10%-20% and may be only a few percent. We analyze the clumping 
property of the hydrogen ionized by QSOs and study the evolution of the Stromgren sphere. Wc find 
that the expected Stromgren radii from our models are consistent with observations if the lifetime of 
the highest redshift QSOs is about or longer than a few times 10 7 yr (as is the lifetime of the main 
population of QSOs; with comoving number density peaked at z ~2-3). With such a QSO lifetime, the 
ages of most of the observed QSOs are long enough that the QSO photon emission is balanced by the 
recombination of the hydrogen ionized by QSO photons in their Stromgren spheres, and the expected 
Stromgren radii from the balance are independent of the detailed values of the QSO ages. We also 
point out a statistical method involving a larger sample of QSOs having Gunn-Peterson troughs in 
future observations, which may potentially check or rule out the possibility that the highest redshift 
QSOs have a shorter lifetime (e.g., < 10 7 yr), even without an accurate estimate of the hydrogen 
clumping property. 

Subject headings: cosmology - theory: cosmology - early universe: galaxies - formation: galaxies: 
high redshift - quasars: general 
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1. INTRODUCTION 

QSOs at redshifts z > 6.1 have been discovered 
by the Sloan Digital Sky Survey (SDSS) recently 
l|Fan et al]l2001l 120031 12004(1 ■ Their spectra have shown 
the existence of the Gunn-Peterson ab sorption trough 
(|Ounn fc Peterson! 119651: iScheuerl 119651) at wavelengths 
blueward of the Lya emission line, which suggests that 
hydrogen in the earl y universe is significantly neutral (> 
1% in mass average: | B e cker et al.ll2001tlDiorgovski et al.l 
120011: iFan et all I2001L 12002^ Despite the presence of 
the Gunn-Peterson trough in the spectra, a significant 
amount of flux is transmitted at wavelengths blueward 
of but near the center of the Lya line, which suggests 
that the regions surrounding these highest redshift QSOs 
are highly ionized by QSO photons (i.e . , the proxim- 
ity effect of QSOs, e.g., in iBaitlik etaD 11988ft . These 
highly ionized regions are generally idealized as spheri- 
cal and called "Stromgren spheres." Their radii (or the 
Stromgren radii) can be estimated from the wavelength 
difference between the Gunn-Peterson trough and the 
Lya line center fe.g. lWhite et alJl2003ft . 
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Physically the sizes of the highly ionized regions are 
closely related to not only the neutral hydrogen distri- 
bution in the local universe of these QSOs before their 
nuclear activities turn on but also the subsequent QSO 
luminosity evolution, either of which thus may be con- 
strained by the Stromgren radii estimated from QSO 
spectra provided that the range of the other one is 
given. For example, using the observationally deter- 
mined Stromgren radii and assuming the neutral hydro- 
gen fraction (defined by the ratio of the neutral hydrogen 
number density to the total hy drogen number density) 
,thi = 1. lHaiman fc Cer] (|2002f) found that the lifetime 
of one highest redshift QSO is not sh orter than 2 x 10 7 yr 
(see also Pentericci et al. 2002) and iWhite et all l)2003ft 
got a substantially shorter lifetime for two; using the 
measured Stromgren radii of two highest redshift QSOs 
and some independent con straints on the QSO life time 
obtained by other studies, IWvithe fe Loebl l(2004aft ar- 
gued that xhi may be higher than tens of percent at 
z ~ 6.3, which is significantly higher than the lower 
limit of cchi ^ 1% constrained by the flu x upper limit 
of the Gunn-Peterson absorption trough l|Becker et alJ 
200x1: IFan et al J 12002) IWhite et al.ll200l In ad dition, 
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tigated the structure of the transmitted spectra at wave- 
lengths blueward of the Lya line and pointed out that the 
red damping wing of the Gunn -Peterson trough might 
leave detectable features there ( Miralda -Escudel ll9981 . 
wh ich may also help to c onstrain xm (see a recent study 
by iMesinger et al J 120041 using mock data from simula- 
tions. 

In most previous investigations on the reionization 
process around the highest redshift QSOs, the QSOs 
are implicitly assumed to be located in cosmic av- 
erage regions. However, in the hierarchical forma- 
tion and evolution scenario o f galaxies and QSOs (e.g., 
iKauffmann fc Haehneltl 120001) . luminous QSOs (espe- 
cially those at high redshifts z > 6.1) are located in 
rare overdense regions. The structure formation in these 
regions and consequently the gas distribution, star for- 
mation, and reionization of the neutral hydrogen may 
be significantly biased away from those in the regions 
with cosmic mean density. Some discussions in a few re- 
cent papers are also rela ted to some biasing effects of a 
dense environment (e.g .. lCenll2003bt iCi ardi et all 120031 
[ Furlarietto et all 120041 iGnedin fc Pradal 120041 ISantosI 
120041 iWvithe fc Loebl l2004bfl . The purpose of this pa- 
per is to quantitatively investigate the dense environ- 
ment of the highest redshift QSOs, the reionization in 
their local universe due to ionizing photons from both 
stars and QSOs, and the evolution of their Stromgren 
spheres. Comparison of the expected Stromgren radii 
from the more realistic model with observations would 
improve our understanding of the universe's reionization 
history and the formation and evolution of QSOs. 

This paper is organized as follows. The observational 
properties of the highest redshift QSOs that have the 
Gunn-Peterson trough are reviewed in § In § we 
investigate the environment around these QSOs in the 
framework of the cosmic d ensity perturbation a nd struc- 
ture formation theory fe.g.. lLacev fc Colell993[) . We first 
estimate how dense the environment is from the number 
density of the highest redsh ift QSOs. Then we int roduce 
a hybrid model (see also iBarkana fc Loebl l2004j) based 
on the excursion s et app roach of the Press-Schechter for- 
mula l)Bond et al.lll991t lLacev fc Colelll993|) to see how 
the structure or halo formation is enhanced in these rare 
overdense regions in § 13.21 Gas in halos may be shock 
heated, collisionally ionized, and then cool to form stars, 
galaxies, and/or QSOs. Ionizing photons from stars and 
QSOs will jointly contribute to the reionization of the 
surrounding neutral hydrogen. In this paper we iso- 
late their effects and first consider the reionization due 
to stars (the variable cchi in this paper represents the 
neutral hydrogen fraction obtained after considering the 
reionization due to stars but before the reionization due 
to the QSO becomes effective). The justification for this 
isolation is that before the cosmic time that the QSO 
redshift corresponds to, the QSO ionizing photon emis- 
sion rate was high enough to be comparable to (or higher 
than) the ionizing photon emission rate from stars only 
in a short period (a few times 10 yr) and the cumula- 
tive number of ionizing photons are mainly from stars 
(see detailed discussion in § 13.31 and 14. 1(1 . In S 13.31 using 
some simple star formation models, we find that more 
than 80%-90% of neutral hydrogen in these dense re- 
gions is ionized by photons from stars. Inhomogeneity 
or dumpiness of hydrogen may significantly affect the 



reionization process. In § 13.41 we give a rough estimate 
of the dumpiness of hydrogen ionized by QSO photons, 
based on the halo distribution obtained in § 13.21 some 
numerical simulation results obtained in other studies, 
and some physical arguments. In the environment of the 
highest redshift QSOs analyzed in §13 we study the evo- 
lution of the apparent size of the Stromgren sphere highly 
ionized by QSO photons in § 0] The model results are 
compared with observations and the implications are dis- 
cussed. The main conclusions are summarized in § [SJ 

Throughout this paper, we adopt a standard ACDM 
cosmological model with 17a = 0.73, Q m = 0.27, fib — 
0.0444, ff n = Jlkm s" 1 Mpc -1 , cr 8 = 0.84, and n = 0.93 
(|Spergel et al.ll2003ft . 

2. SAMPLE OF THE HIGHEST REDSHIFT QSOS THAT 
HAVE GUNN-PETERSON TROUGHS 

To date, four QSOs have been observed at redshift 
z > 6.1, and all of them are detect ed to have the 
Gunn-Peterson trough in their spectra l)Fan et alJl2003t 
120041: iWhite et all l2003() and are compiled in Tabled 
As seen from Table the redshifts of the observed 
QSOs are in the range of zq ~6.2-6.4. The redshift 
may be dete rmined from the Lya line (e.g., for SDSS 
J1623+3112: iFan et all 120041) . and the typical redshift 
error is ±0.02, which is mainly caused by the uncertainty 
in the determination of the location of the line center in 
the spectra. The redshift of SDSS J1148 +5251 is pre- 
cisely determined from C O and Mgn lines l)Walt,er et alJ 
120031 IWillott et aT]l2003D . and we simply assign a small 
error of ±0.005 to its redshift (similarly for SDSS 
J1048±4637 whose r edshift is determined by Mgu line; 
iMaiolino efaTl 120041) . Although the redshift of SDSS 
J1030+0524 has also been determined from metal lines 
(e.g., Civ), the redshift error can still be as large as 0.02, 
caused by the possi ble blueshift of the Civ line (e.g., 
iRichards et al.l2002[) , and we therefore set its redshift er- 
ror to be +0.02 and —0.005. The Stromgren radii of three 
QSOs in Tabic ^ have been estimated through the dif- 
ference between the QSO redshift and the redshift of the 
onset of the Lya Gunn-Peterson trough in their spectra 
in the literature (i.e., knobs' see a ^ so z mi,obs in Tab.^ 
for the corresponding redshift range of the Stromgren 
sphere, except for SDSS J1048±4637). The errors on 
r Hii obs i n Tabic H a rc introduced by the uncertainty in 
the determination of these QSO redshifts. 

There may also exist uncertainties in the determina- 
tion of the onset redshift of the Gunn-Peterson troughs, 
which further introduce uncertainti es in the estimate of 
the St romgren radii. For example, IMesinger fc Haimanl 
(2004) point out that the radius estimated through the 
Lya trough may only give a lower limit of the Stromgren 
radius and that a better estimate may be obtained 
by using the Ly/3 trough. Using the Ly/3 trough and 
some numerical simulation results, IMesinger fc Haimanl 
(2004) increase the estimate of the Stromgren radius of 
SDSS J1030±0524 to - 6Mpc (from - 4.6 Mpc listed 
in Tab. 1). Note that the new estimate might be af- 
fected by the numerical sim ulation used in the analysis 
(|Mesinger fc Haimanl 12004)) . which does not well repre- 
sent the overdense environment of the highest redshift 
QSOs discussed in this paper. In addition, other QSOs 
listed in Table ^ do not have the Stromgren radius es- 
timated by using the Ly/3 trough in the literature, and 
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thus we still mainly use the values estimated from the 
Lya troughs for consistency in this paper and discuss 
the implication of a possible systematically larger value 
of the observational Stromgren radius in § 14.31 

We estimate the ionizing photon emission rates of these 
QSOs, Np hs q, by using their absolute magnitude M1450 

at 1450 A listed in Table ^ (compiled from iFan et all 
1200 1L 120031 12004) and the spectral templ ates of rela- 
tively lo w redshift QSOs {z < 3) shown in iTelfer et alJ 
(2002). ijFan et alJ [20041 have found that the average 
sp ectrum of the high est r edshift QSOs is sim ilar to that 
in ITelfer et ail 120021 ) In ITelfer et all l|2002fl . the aver- 
age spectra of the low-redshift QSOs (energy radiated 
per unit time and per logarithmic interval of frequency 
vf v oc v 1+a ) is well fitted by a double power law in the 
optical to EUV band. For radio-loud QSOs, the EUV 
and near-UV spectral indices are oeuv = —1-96 ± 0.12 
and «nuv = — 0.67±0.08, respectively, with the break at 
~ 1280A, while for radio-quiet QSOs they are oeuv = 
-1.57 ±0.17 and a NUV = -0.72 ±0.09, respectively. We 
count photons with energy in the range 13.6-54.4 eV as 
ionizing photons. Since so far little information on the 
radio properties of the highest redshift QSOs is avail- 
able, we obtain two emission rates for each QSO listed 
in Tabled the smaller of which is obtained by adopting 
the spectral indices of radio-loud QSOs and the larger 
of which is obtained by adopting the spectral indices 
of radio-quiet QSOs. We take the middle point of the 
whole range of the two values (including their errors) as 
iVp hs q and take the whole range as the error range of 
iVp hs q (see also the filled circles and their horizontal er- 
ror bars in Fig. Et). Note that for SDSS J1030±0524, 
compared to the rate listed in Table HI a substantially 
smaller rate is obtained in Mc singer fe Haimanl ( 2004) 
by fitting the spectrum of the QSO in a surrounding hy- 
drogen density distribution given by a numerical simu- 
lation. As mentioned above , the n umerical simulation 
used in iMesinger fc Haimanl <)2004[) mimics the cosmic 
average environment rather than the rare overdense en- 
vironment around the highest redshift QSOs; thus the 
smaller rate may be an underestimate since more ionizing 
photons would be required to balance the recombination 
in a denser environment. 

The highest redshift QSOs are extremely rare and their 
comoving number density Nq is roughly the same as the 
density of the QSOs with a bsolute magnitude A£i 450 < 
-26.7 at z ~ 6 estimated bv lFan et all (|2003l 127)041) . i.e., 
N Q ~ (6± 2) x 10 -10 Mpc" 3 . (Note that the value of N Q 
is obtained by assuming a cosmological model with Hq = 
65kms _1 Mpc -1 , Q m = 0.35, and Qa = 0.65; we have 
adjusted this value to the cosmological model adopted in 
this paper and found no significant difference.) 

3. THE ENVIRONMENT AROUND THE HIGHEST 
REDSHIFT QSOS 

The highest redshift QSOs are located in rare over- 
dense regions according to current structure and QSO 
form ation and evolution picture in t he ACDM cosmogony 
('e.g-.TKauffma nn fc Hae hnclt 200fj). To quantify the en- 
vironment of these QSOs, we first estimate the mean 
mass overdensity in their surrounding regions in § 13.11 
Then we study the structure and halo formation in these 
overdense regions in § 13.21 and study the star forma- 



tion and the reionization due to ionizing photons emitted 
from these stars in 5 13.31 In § 13. 41 we study the dumpiness 
of hydrogen in the Stromgren sphere of the QSOs. 

3.1. Overdensity 

The mean overdensity within a_sphere with proper ra- 
dius r at redshift z is defined by 5 r = \(p) r — p]/p, where 
(p) r is the average comoving mass density in the sphere 
and p is the comoving cosmic mean mass density. We de- 
note the comoving radius of the sphere by R = r(l + z). 
(Here we do not distinguish the difference between the 
Eulerian radius and the Lagrangian radius, since the 
overdensity considered in this paper is small.) The vari- 
ance of the density fluctuations at redshift z in spheres 
with comoving radi us R, v 2 (R, g ) ; can be determined by 
(e.g., see eq. 2.55 in lWhitdll994|) 

r°° dk 

a 2 (R,z) = D 2 (z) k 2 P(k)W 2 (kR), (1) 

Jo 27r 

where D(z) is the linear growth factor of perturbations 
at z and normalized to 1 at the present time, P(k) is 
the power spectrum of the density pert urbation field and 
is com puted with the fitting formula of lEisenstein fc H7J 
GHU, and W(kR) = 3[sin(fci?) - kRcos(kR)]/(kR) 3 is 
the top-hat window function. If the perturbations on 
the scale of r are still in the linear growth regime at 
redshift z, the mean overdensity S r follows a Gaussian 
distribution: 



p(5 r , z)dS r 
1 



27Rj(i?, z) 



exp 



2a 2 (R,z) 



dS r , 



(2) 



and the probability that a random region with scale r 
has a mean overdensity higher than a value 5' r is given 

by 



V{> S' r ,z) = J p{S r ,z)dS r 



(3) 



Since the Stromgren radii measured from the observed 
QSO spectra are around 5 Mpc, below we consider a 
sphere with proper radius rQ — 5 Mpc. According to 
equation QJ, the standard deviation of the density fluc- 
tuations in this sphere at z ~ 6 is a[R — ?*q(1 + z), z ~ 
6] ~ 0.06, which is much smaller than the critical over- 
density 8 C = 1.68 for dynamical collapse, and the per- 
turbations are still in the linear growth regime. Thus 
equations (0) and Q can be applied to these spheres. 

In regions with sufficiently high overdensity, suffi- 
ciently massive halos will form and these halos are the 
hosts of luminous QSOs. In a random spherical region 
with proper radius tq = 5 Mpc at z ~ 6, the probability 
of finding a halo capable of hosting nuclear-active QSOs 
with Mi4 50 < -26.7 (e.g., with mass ~ 10 13 M Q ) is 



Am- 



4 t-1 
duty 



(4) 



where /duty (z) represents the duty cycle of QSOs, defined 
to be the number ratio of the nuclear active QSOs to the 
halos capable of hosting the QSOs at redshift z. The 
halos capable of hosting QSOs include both those host- 
ing a nuclear-active QSO and those hosting a dead QSO 
with nuclear activity quenched. For the highest redshift 
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QSOs, we simply assume /duty ~ 1 because most suffi- 
ciently massive halos (e.g., > 10 13 M Q ) at z ~ 6 are ac- 
tually fo rmed within a peri od (from z = 7 to 6; e.g., see 
Fig. 1 in lMo fe Whitdl2002|) not significa ntly longe r than 
the QSO lifetime (> a few times 10 7 vr: lYufcLull2004al: 
lYu fc Tremaind l2002t see more discussions about the 
QSO lifetime in §0J. By multiplying equation © by a 
factor of 2 (this factor comes from the so-called clou d-in- 
cloud problem; see discussions in Bond ct al. 1991) and 
then setting it equal to equation we obtain that the 
mean overdensity of the regions associated with the high- 
est redshift QSOs is 5 rQ > S' rQ ~ 4a (R, z ~ 6) ~ 0.25. 

Since the overdensity distribution g(S r ) decreases expo- 
nentially at the high overdensity end, most of these over- 
dense regions associated with the highest redshift QSOs 
do not have a mean overdensity significantly higher than 
S' rQ , but near the value. 

Note that the connection of equation © with equa- 
tion above has assumed a one-to-two correspondence 
between the likelihood of finding a region with mean 
overdensity above a certain value and the likelihood of a 
same-sized region containing a halo capable of hosting a 
quasar. Although in principle there might exist regions 
with overdensity above the certain value but not con- 
taining such a halo, the simple approach above is good 
enough in practice for the rare overdense regions and for 
the purpose of this paper. For example, given a halo with 
a certain mass, the method in Barkana (2004; for sim- 
plicity the details are not presented here), which is also 
based on the (extended) Press-Schechter formalism, may 
provide the expected average overdensity within a certain 
region around the halo and the standard deviation of the 
overdensity. Assuming that the host halos of the high- 
est redshift QSOs are aro und 10 13 M W (o r 2 x 10 13 M Q ) 
and using the method in Barkana (2004), we obtain an 
average overdensity of ~ 0.13 (or 0.16) within the region 
with radius tq surrounding the halo, which is smaller 
than the estimate above. Hereafter, we simply adopt 
& 



J r Q 



S' rQ ~ 0.25 as the mean overdensity of the spheres 



with proper radius rq> = 5 Mpc centering on the highest 
redshift QSOs. Using the value of 5 rQ = 0.13 will not 
qualitatively affect the conclusions in this paper. 

3.2. The formed structure 

In this section, we introduce a hybrid model to esti- 
mate the mass distribution of halos formed in the rare 
overdense regi ons associated with th e highest redshift 
QSOs (see also lBarkana fc Loebll200l . 

We define the halo mass function as jjj(z) so that 
^jjdM represents the comoving number density of halos 
with mass in the range M — > M + dM at redshift z. In 
the cosmic average regions, the halo mass function can 
be given by 



dN 
dM 



P_ 

M 



da 2 (R c , z) 



dM 



f[S c ,a 2 (R c ,z)}, 



(5) 



where R c = (SM/Airp) 1 / 3 , and f[6 c , a 2 (R c , z)]da 2 (R c , z) 
is the mass fraction of halos with mass in the range M — > 
M + d M at redshift z. In the model of lPress fc Schechterl 
ltl97l - f[5 c ,a 2 (R c ,z)} is given by 

1 v . i > 

/2tt a 2 (R c ,z) 2 



fpa[6 c , o 2 (R c ,z)} 



where v = 5 c /o~(R c , z). Compared to the results obtained 
from cosmological simulations with very large volumes, 
however, the Press-Schechter mass function substantially 
underestimates the abundance of the rare massive halos 
that host galaxies at high redshifts and overestimates 
the abundance of intermediate-mass h alos. To more ac- 
curate ly match the simulation results, iSheth fc Tormenl 
(1999) introduce a new formula 



f ST [S c ,a 2 (R c ,z)] 

A' Vo 1 ^ 
2tt a 2 (Rc,z) 



{a'v 2 )i 



exp 



a'v 2 



(7) 



with the best-fit parameters a! = 0.707 and q' = 0.3, 
and the normalization parameter A' = 0.322. This for- 
mula can be derived by introducing an ellipsoidal col- 
lapse model instead of the usually adopted spherical col- 
lapse model into the excursion set appr oach of the Press- 
Sc hechter formula (IShet h et alJ f2001j) . In the model 
of ISheth fc Tormenl 119991. we have the cosmic average 
mass fraction of halos with mass higher than M given 

by 



i ? ST,halo(> M,Z) = 



f ST (S c ,a 2 )da 2 



(8) 



In a rare overdense region with volume much smaller 
than the universe, the halo mass distribution significantly 
deviates from the cosmic mean distribution shown by 
equations I©-©. Similar to using the excursion set 
approach to get t he extended Press- Schechter formula 
l|Bond et alJll99lHLacev fc Coldll993ft . we have the cor- 
responding biased mass fraction of fps[S c , c 2 (Rc, z)] in a 
region with comoving radius R and mean overdensity 5 r 
at redshift z given by 



f^r[S c ,a 2 (R c ,z);6 r ,R] 
= fps[S c - S r ,a 2 {R c ,z) - 



a 2 (R,z)}, 



(9) 



and the mass fraction of halos with mass higher than M 
given by 



F^Uo(>M,z) = 



erfc 



^2^(R c ,z)-a 2 (R,z)} 



(10) 



As shown in iBarkana fc Loebl l|2004D . the halo mass 
function in overdense regions can be obtained by ad- 
justing the Sheth-Tormen formula with a relative cor- 
rection bas ed on the extended Press-Schechter model 



l|Bond et alll99lHLacev fc Colelll993D . which can match 
cosmological simulation results well. Similar to this hy- 
brid model, the biased halo mass function in the regions 
with mean overdensity S r can be obtained by 2 



dN' 

dM_ 
da 2 (R c ,z) 



dM 



(1 + S r )p 
M 

f h ™{8 c ,a 2 (R Cl z);5 r: R} : 



(11) 



2 The biased halo mass function in a high-density region may 
also be obtained by using the e llipsoidal collapse moving barrier 
model in Shcth fc Tormenl 120021) instead of the hybrid model used 
here. 
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where 

/ bias [<5 c , a 2 (R c , z); S r> R] = f ST [S c , a 2 (R c , z)\ 

f^ s [Sc-~S r ,<T 2 (Rc,z)-a 2 (R,z)} 

fp S [S c ,a 2 (R c ,z)} ■ [ } 

In these overdense regions, the mass fraction of halos 
with mass higher than M is given by 

F^(>M,z) = 

f hi ^(6 c: a 2 ;S r ,R)da 2 . 



(13) 



In halos with sufficiently high masses, most of the gas 
in them is collisionally ionized and some of the ionized 
gas may cool to form stars, galaxies, and/or QSOs. Neu- 
tral hydrogen in the universe may then be photoionized 
by photons from stars and QSOs. We define a character- 
istic mass scale M m ; n , the minimum mass of halos whose 
potentials are deep enough to retain ionized gas in the 
photoionization equilibrium. The temperature of ionized 
gas in photoionization equilibrium is ~ 10 4 K, and we 
have M min = 3.4 x 10 7 M Q [10/(1 + z)} 1 - 5 (see details for 
the determination of this mass scale in Appendix B in 
iBenson et all 120071) . We include both the gas outside 
halos with mass greater than M m j n and that inside these 
halos but not collisionally ionized as IGM (the mass of 
baryons in galaxies is negligible; see Fig. ^ below). In 
these regions, the mass fraction of hydrogen that is in 
IGM is given by 



xh,igm(z) 



M 



(1 + S r )p 



dN 

dM 



-i bias 



dM, (14) 



where yn is the fraction of hydrogen that is collisionally 
ionized in a halo with mass M at the virial tempera- 
ture Tvi,., and the ratio of the baryonic mass to the total 
mass in halos is assumed to be the same as the cosmic 
mean fl h /Q m . The virial temperature is T vir = ^f^f^, 
where /i ~ 0.6 is the mean molecular weight, mn is the 
proton mass, ks is the Boltzmann con stant, and r v \ r is 
the v irial radius of the halo (see, e.g.. iBarkana fc Loebl 
|2001|) . The determination of the detailed distribution 
of collisionally ionized hydrogen (see also eq. EHI below) 
would need to take into account detailed gas dynamics 
inside the halos (e.g., photo-evaporation processes, dy- 
namics related to galactic disks), which is currently not 
well understood. The t/h us ed in this paper comes from 
iSutherland k, Dopital l|1993t) , as commonly adopted. 

Our calculation shows that in the rare overdense 
regions associated with the highest redshift QSOs, 
£h,igm (tq , zq ~ 6.2 — 6.4) ~ 0.88 — 0.89 while the cosmic 
mean xh,igm(zq) is around 0.94. 

3.3. The reionization due to stars 

Star formation in the formed halos is directly related 
to the cooling of their virialized gas due to radiation by 
atomic, molecular, or metal lines. Gas cooling due to 
atomic lines, and thus star formation, can be very ef- 
ficien t in halos with T v , r > 10 4 K, as usually argued 
(e.g.. IHernauist fc Springeil f2003). However, atomic line 
cooling is not efficient in halos with T vir < 10 4 K (the 
so-called minihalos). Molecular hydrogen (H2) line cool- 
ing may be important for the very first (Population III) 



star formation in those minihalos at very high r edshift 
(e.g., z > 10; see lAbel et aTll2002t iBromm et al.lll999jh 
These Population III stars may lead to the reionization 
of the universe at very high redshift, but this may be 
unimp ortant for th e (second) reionization that finished at 
z ~ 6 (Cen 20 03al) . At lower redshift, external UV radia- 
tion may suppress the molecular H2 abundance in those 
minih alos and thus redu ce the cooling due to H2 lines 
(e.g.. lHaiman et aTlll996t) . However, the precise cooling 
threshold and H2 abundance in these regimes are not well 
known. Furthermore, star formation in minihalos is also 
likely to be suppressed by the processes of photoevapora- 
tion and supernova disruption. Below we only consider 
star formation in halos with T v ; r > 10 4 K and neglect 
star formation in minihalos. Detailed physical processes 
of star formation in halos with T vu > 10 4 K are com- 
plicated and not fully understood. Here we adopt two 
simple models to estimate how many baryons form stars 
in these halos: 

• Model (a): The star formation efficiency (i.e., the 
mass fraction of baryons to form stars) in halos 
is assumed to be i n depend ent of the halo mass, 
as adopted in iCenl (|2p03a) . To satisfy the con- 
straint that the reionization epoch ends at z ~ 6, 
ICenl l)2003al) obtained a best-fit value of 0.1 for 
the star formation efficiency in halos with a virial 
temperature greater than 10 4 K. We adopt this 
value and thus have the mass fraction of baryons 
that are in formed stars in the rare overdense re- 
gions associated with the highest redshift QSOs, 
g^(z) ~ 0.1F l ;X(> Mi, z), where M 4 is the mass 
of halos with virial temperature of 10 4 K. 

• Model (b): the star formation efficiency is assumed 
to depend on the halo mass. After considering the 
feedback due to galactic winds in the sta r forma- 
tion process, [Hcrnquist & Springe] (2003) argued 
that the normalized star formation rate (the mass 
fraction of baryons to form stars per unit time) 
may maintain a roughly constant level for halos 
with virial temperatures in the range 10 4 -10 6 ' 5 K 
but may rise to a level approximately 3 times higher 
for halos with hotter virial temperatures. Based on 
this approximation, the star formation rate in ha- 
los with temperature T v -, r is given by (see eqs. 41 , 
43 and more details in lHernauist fc Springell l2003') 

S(T V i I} z) ~ 



s q(z) 

3s q(z) 





for 10 4 K < T vir < 10 6 5 K 



for T vir > 10 
otherwise, 



6.5 



K. 



(15) 



where q(z) 
«a] 1/3 



9/27) 



x = [n m (i+z) 3 + 



m = 6, and rj = 1.65. For the cosmological 
model adopted in this paper, we adjust \ from 4.6 
to 4.4 and s from O.OO6/1 Gyr" 1 to O.OO8/1 Gyr" 1 
(see details for the depend ence of y and sg on 
cosm ological parameters in iHernauist fc Springe! 
2003) . Thus the mass fraction of baryons that are 
in formed stars in the rare overdense regions asso- 
ciated with the highest redshift QSOs is given by 



s q(z) x 



6 



{F^? (> M 4 , z) + 2F h b X(> Ma. 6 , *)]-<&, (16) 
while the corresponding cosmic mean is given by 

/>OC 

ffstar^) — / S q(z) X 
J z 

[FsT,halo(> M 4 ,Z) + 2F ST ,halo(> M 6 . 5 , z)]^-dz ,(17) 

where M6.5 is the mass of a halo with virial tem- 
perature 10 6 - 5 K. 

We plot the results of g^l^(z) and g s tar(z) in Figure ^ 
for both models a and b. Their values at the QSO red- 
shift zq ~ 6.2 — 6.4 are summarized in Tabled As seen 
from Figure 12 both models indicate that star formation 
in the overdense regions (solid line) is significantly bi- 
ased away from or higher than the cosmic average at 
the same redshift (dashed line). Although the mass frac- 
tions of baryons in formed stars obtained from models a 
and b are not exactly the same, the two models give a 
consistent result in that the mass fraction of baryons in 
formed stars in the regions associated with the highest 
redshift QSOs at zq ~ 6.2 — 6.4 is similar to the cor- 
responding cosmic mean values at a later time z ~ 5, 
i.e., 5star( z Q) — 9star(z ~ 5) (indicated by the arrows in 
Fig. ^ see also Tab. [2J. Since the number of ionizing 
photons emitted from stars is proportional to the mass 
of baryons in formed stars, the cumulative number of the 
ionizing photons produced in the stars in the rare over- 
dense regions associated with the highest redshift QSOs 
at zq ~ 6.2 — 6.4 is thus similar to that of the correspond- 
ing cosmic mean at a later time z ~ 5. However, their 
hydrogen ionization statuses may not be naively simi- 
lar, for example, because of differences in their proper 
densities or in the process/time for hydrogen to dynami- 
cally respond to the ionization (especially for the photo- 
evaporation process of minihalos). 

We assume the initial mass function for star forma- 
tion in the surrounding regions of the highest redshift 
QSOs is sim ilar to that measured in the nearby universe 
l)Scalolll998tl and thus that roughly 4000 ionizing pho- 
tons are produced per baryon in the formed stars (e.g., 
iBarkana fc Loebll2001|) . The rate of ejection of the ion- 
izing photons into the IGM can be given by 

Att dn hlas (7\ 
Afeax * Y R3{1 + ^ JS ^ Ll X 400 °/esc>(18) 

where iVb is the baryon comoving number density and 
/ csc is the fraction of the ionizing photons escaping into 
the IG M and is roughly in the range 0.1-0.2 (e.e.. ICenl 
l2003a|L In spheres with proper radius tq = 5Mpc 
associated with the highest redshift QSOs, we obtain 
-ZV^h-fstar^Q) for both models a and b and list them in 
Table [21 For comparison, the cosmic mean injection 
rate of ionizing photons due to stars at z ~ 5 in the 

same comoving volume, iVphs.star — ^-R 3 Nb dSat ^ z ^ x 
4000/csc, is also listed in Table [3 which is roughly 
half of Np^ staT (zQ). Assuming that the background 
ionizing energy intensity 3 V (v. photon frequency) is 
mainly contributed by stars, we estimate that an ion- 
izing photon emission rate of -/Vphs.star leads to the 
background intensity J -21 at the Lyman limit in units 



of 10" 21 ergs cm" 2 s" 1 Hz" 1 sr" 1 of about 0.06-0.13 
for model a and .1-0.2 for model b (using eq. 29 in 
iMadau et alj UggO 1 ). These numbers are roughly consis- 
tent with tho se determined fr om observations at z ~ 5 
(see Fig. 2 in lFan et al.ll2002T) . which suggests that the 
star formation models used in this paper are roughly con- 
sistent with observations. In the rare overdense regions 
around the highest redshift QSOs, according to equation 
1)18(1 the cumulative number of the ionizing photons emit- 
ted from stars before t(zQ) is about 5-10 per hydrogen 
atom in the IGM within the Stromgren sphere, which 
is substantially higher than that from the central QSO 
(< 1 — 2; see also discussion in S I4.1JI . 

As shown above, a large number of stars have already 
formed in the rare high-density regions before z ~ 6, and 
in most of the time before the cosmic time t(zQ), stars 
are the main sources of the ionizing photons contribut- 
ing to the re-ionization of neutral hydrogen in the IGM 
(as mentioned in the introduction and as will also be dis- 
cussed in § 14.1(1 . The reionization starts in low-density 
regions (or voids) and then gradually pen etrates deeper 
into high-density regions (se e discussions in lGnedinl2000t 
iMiralda-Escude et al.l l20001 . Thus, before the reioniza- 
tion by the QSO photons becomes effective, the IGM in 
low-density regions surrounding those stars (or galactic 
units) has already been highly ionized, while the IGM 
in high-density regions may still be significantly neutral 
since the ionizing photons from stars are still insufficient 
to ionize all the IGM. The division from the highly ion- 
ized regions to the significantly neutral regions may be 
characterized by a critical overdensity A cr i t , where the 
overdensity A is defined through A = pb/pb, Pb is the 
comoving baryon mass density at a given space posi- 
tion, and pb is the cosmic comoving mean mass den- 
sity of baryons (note the difference from the definition 
of the overdensity 5 r , which is averaged over a proper 
scale of r). After the nuclear activity of the QSO turns 
on and the emission of its ionizing photons becomes sig- 
nificant, QSO photons will ionize some remaining neu- 
tral hydrogen (in high-density regions) not ionized by 
stars. Compared to ionizing photons from stars, the pho- 
tons from QSOs are much harder and can penetrate into 
high-density regions more deeply. In addition, QSOs are 
rare point sources located in high-density regions, while 
stars are embedded in relatively numerous halos and act 
like diffuse sources. A significant fraction of the ioniz- 
ing photons from stars may leak out of the regions (with 
tq = 5 Mpc) considered in this paper, while ionizing pho- 
tons from QSOs do not, but are absorbed. 

Here we give a rough estimate of the ionization status 
due to stars in the overdense regions associated with the 
highest redshift QSOs, such as the critical value A cr it 
and the distribution of A: 

• We assume that the IGM is in the local photo- 
ionization-recombination equilibrium, as com- 
monly adopted (see, e.g.. lFan et al]l2*002jk that is, 
the photoionization rate nuiT is balanced by the 
recombination rate nmin e aB(T), where nni, Uhii, 
and n e are the proper number densities of neu- 
tral hydrogen, ionized hydrogen and electrons in 
the IGM, respectively, T is the number of pho- 
toionizations per H atom per unit time and is 
roughly proportional to the ionizing photon emis- 
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sion rate, and ojb(P) is the recombination coeffi- 
cient at temperature T l|Abel et al.lll997T) . Note 

that A^ star (zq ~ 6.2 - 6.4) ~ 2AV phS;Star (z ~ 5) 
and that the mean hydrogen density in the rare 
overdense regions associated with the highest red- 
shift QSOs at zq is higher than the cosmic mean 
at z ~ 5 by the same factor of 2. Thus, compared 
to the photo- ionization equilibrium at z ~ 5, given 
the same ionization ratio of riHii/^Hi, the effect of a 
higher ionizing photon emission rate at z ~ 6.2—6.4 
is compensated by the effect of a higher hydrogen 
density in the rare overdense regions. Hence, we set 
the critical overdensity A cr ;t within the rare over- 
dense regions associated with the highest redshift 
QSOs at zq ~ 6.2 — 6.4 to be similar to that in the 
cosmic average regions at z ~ 5. The cosmic aver- 
age ionization state of the IGM a t redsh ift z ~ 5 
has been measured by iFan et all (|2002fl through 
the high-redshift QSOs around that redshift an d 
Acrit — 10 — 20 (see Fig. 8 in IFan et alJ 120021. 
The simple arguments above have circumvented 
a more sophisticated and lengthy way to obtain 
Acrft) i.e., using N^ Ktar (z Q ) and repeating some 
detail ed formalisms used in iMiralda-Escude et al.l 
(200$, which gives an even larger A cr i t . Hereafter, 
we simply adopt A cr j t ~ 10 — 20. At least the 
regions with overdensity below this value can be 
highly ionized by stars. This value is also compati- 
ble with the constraint on the rcionization process 
from the Stromgren sphere in § 0J 

• Since the mass fraction of halos with mass greater 
than M m ; n within the rare overdense regions asso- 
ciated with the highest redshift QSOs is similar to 
the cosmic mean fraction at a later time z ~ 4.5 — 
4.8 [see Fig. UK, where F halo (> M 4 , z) ~ 10g stax (z), 
or see Tab.0, we assume that the volume- weighted 
mass density distribution of baryons outside of ha- 
los in these rare overdense regions is analogous to 
(or at least does not significantly deviate from) 
the cosmic mean distribution at a later time z ~ 
4.5 — 4.8. Numerical simulations have shown that 
the cosmic mean volume- weighted density distribu- 
tion function P(A), defined so that P(A)dA rep- 
resents the fraction of the volume occupied by re- 
gions with over density in the range A — > A + dA, 
is described by ijMiralda-Escude et alJl^O OO) 



P(A) = Aexp 



(A" 



-2/3 



Co) 



2 1 



2(V3) 2 



A _/3 , 



(19) 



where A and C'o are set by normalizing 
J P(A)dA = 1 and / AP(A)dA = 1, and 8 and [3 
are given by simulations, for example, So — 3 and 
(3 = 2. 41 at z ~ 5 and 6 = 3.6 and (3 = 2.25 at 
z ~ 4 (jChiu et al.l[2 003'). The P(A) given by simu- 
lations may be biased away from the reality because 
of the small size of the box adopte d in numerical 
simulations (|Barkana fc Loebl 12004). but the bias 
is not much for the high redshift zq ~ 6.2 — 6.4 
considered in this paper. 

Using the distribution P(A) given by numerical simu- 
lations at both z ~ 4 and z ~ 5 (eq. I19fl . we obtain 



that the mass fraction of baryons or hydrogen in the re- 
gions with overdensity A < A cr ;t ~ 10 — 20 is about 

x H (A < A^t) = Jo" 1 "" 20 AP(A)dA ~ 0.7 - 0.8. Con- 
sidering that the mass fraction of hydrogen in IGM is 
~ 0.9 in the regions associated with the highest redshift 
QSOs as obtained in § 13.21 and that not only the IGM 
with low densities (A < A cr it ~ 10 — 20) but also part of 
the IGM with high density (A > A cr ;t) may also be ion- 
ized by photons from stars, the neutral hydrogen fraction 
is x m ^ ihjgm - ^h(A < A cr it) < 0.1 - 0.2 in those 
rare overdense regions before the reionization due to the 
QSO becomes effective. 

Not only the models of structure and star forma- 
tion above suggest that a significant number of stars 
have formed in the early universe to contribute to the 
reionization, there also exists emerging observational ev- 
idence suggesting rapid star formation [with a rate of 
~ 3000 M Q yr -1 , which corresponds to an ionizing pho- 
ton emission rate of (0.5 — 0.9) x 10 56 s _1 if the pho- 
ton escaping fraction / csc ~ 0.1 — 0.2] in the host galax- 
ies of the highest redshift QSOs. For example, (1) the 
high infrared luminosity observed in one of the high- 
est redshift QSOs, SDSS J1148+5251, s uggests that a 
large amount of dust has already formed ijBertoldi et alJ 
2003), and (2) the optical study of SDSS J1030+0524 
suggests the existence of super-solar metalic ities in the 
highest redshift QSOs ijPentericci et al. 2002), which are 
also expected to b e the result of rapid star formation 
(Walt er et alj r2003). The very high star formation rate 
suggests that large elliptical galaxies or the host galax- 
ies of those QSOs with masses of 10 11 — 10 12 M Q formed 
on a dynamical timescale of a few times 10 8 yr, which is 
roughly consistent with current scenario of the simulta- 
neous formation of both QSOs/massive black holes and 
elliptical galaxies. If the mass of a massive black hole in 
those QSOs is arou nd a few 10 9 M Q , as estimated by 
iWillott et al. (2003J, the mass ratio of black holes to 
their host galaxies is ~ 10~ 2 — 10~ 3 , which is close to 
the ratio found in nearby massive black holes and galax- 
ies (e .g., iKormendv fc Gebhardd l200lt iMagorrian et alJ 
1998). 

3.4. The clumping factor of the hydrogen ionized by 
QSO photons 

As discussed above, inhomogeneity or dumpiness of 
hydrogen may have a significant effect on the reion- 
ization process. We define a clumping factor Chii = 
((^hii,q) 2 )/(^h) 2 to describe the small-scale dumpiness 
of the hydrogen ionized by QSO photons (not collision- 
ally ionized in halos or ionized by stars), where whii,q 
is the proper number density of the hydrogen ionized 
by QSO photons, and "(• • •)" represents the average 
over the volume of the considered region. Note that our 
definition of the clumping factor here is a little differ- 
ent from that conventionally defined as ((nn) 2 ) / (nn) 2 , 
since the effects of stars and QSOs on the reionization 
are isolated in this paper (see also eq. 1261) . In the re- 
gions with proper radius thii and overdensity 6 rau at 
redshift z, we have (nn)(z) = (1 + <5 rH n)(l + z ) 3 Nh and 
A^h — (1 — 3lHc/4)pb/wH, where helium is assumed 
to be singly ionized and its abundance Yho = 0.24. 
To estimate Chii, we divide the IGM in these regions 
into three components and then consider the contribu- 
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Fig. 1. — Mass fraction of baryons that are in formed stars as 
a function of redshift. The solid line represents the fraction for 
spheres with proper radius tq = 5 Mpc centering on the highest 
redshift QSOs g'' ; "?(2): eg. 1161 . while t he d ashed line represents 
the cosmic mean fraction [§ s tar(z); eq. 1171 . Panels (a) and (b) 
show the results obtain ed f rom different star formation models (a) 
and (b) (described in § 13.31 . respectively. As seen from this figure, 
although the mass fractions obtained from the two models are not 
exactly the same, (zq ~ 6.2 — 6.4) approaches the cosmic mean 
value cjstar(z) of a later tim e z y 5 in both models (indicated by 
the arrows). See details in § 13.31 



tion by each one: (i) Gas in subregions with overdensity 
A < A crit ~ 10 — 20, which has already been fully ion- 
ized by photons emitted from stars, as argued in § 13.31 
This component is not responsible for the absorption of 
the QSO ionizing photons, and thus does not contribute 
to the clumping factor defined above, (ii) The remaining 
neutral gas in halos with mass greater than M m - ln (see 
the definition of M m i n in § I3.2|) , which has not yet been 
collisionally ionized and cannot be photoevaporated out 
of the halo gravitational potential well, (iii) Gas in the 
subregions with overdensity A > 10 — 20 but outside of 
halos with mass greater than M m i n , which has not yet 
been fully ionized by photons from stars. Below we esti- 
mate the contribution by components ii and iii. 

For component ii, it is plausible to assume that the 
remaining neutral gas only in the outer layer of these 
halos is ionized by the photons emitted from stars, and 
most of the remaining neutral gas in component ii will 
be ionized by QSO photons. The average contribution of 
component ii to Chii can be estimated by 



HII 



/int A v 



(i + 4„„) s 



M mln P 



dN 



dM 



bias 



(20) 



where A vir ~ 178 is the virial overdensity and f- mi is 
a weight parameter for the distribution of gas density 
in halos and /j n t — 3.14 if the halo has a nonsingular 
isothermal density profile with a c ore radius of 0.1r v i r 
(see details in § 4 and Appendix B in lBenson et al.l2 001). 
For spheres with proper radius tq — 5 Mpc centering on 
the highest redshift QSOs at zq ~ 6.2 — 6.4, we have the 
mean overdensity 5 rml = 5 rQ — 0.25 (see § I3.1|l . Using 
equations 1)11(1 and 1(20(1. we find that 



The major contributors to C HII are halos with mass in 
the range M m i n -3M mnl and correspond to at most a mass 
fraction of 

1 



(1 + <W) 



(1-Wh) 



M 



dN 



-i bias 



dM 



dM ~ 2%, 



(22) 



of the total baryon. Our calculation also shows that 
the number of halos for ionizing photons emitted from 
the central QSO at a random direction to pass through 
within the scale of tq = 5 Mpc is on average one, which 
is not too small to apply the average clumping factor 
at a random line of sight obtained by equation 1(20(1 to 
the particular observed QSOs. However, for QSOs with 
lower luminosity (e.g., with ionizing photon emission rate 
of less than 10 56 s _1 at z Q ~ 6.2 - 6.4), we find that the 
ionizing photons from the QSO at a significant fraction 
of its lines of sight may not encounter any halo with mass 
greater than M m i n within the Stromgren sphere (but the 
photons could be absorbed by component iii), and thus 
the value of C HII will be an overestimate for the particu- 
lar observed QSOs if they are just located on these lines 
of sight (see discussions in § 14.3(1 . 

Compared to component ii, the contribution of compo- 
nent iii to Chii is not easy to obtain accurately since the 
distribution of the IGM outside of halos and the effects 
of the photoevaporation process on minihalos within the 
rare overdense regions around the highest redshift QSOs 
are less well understood. To get the contribution of 
component iii, we consider the regions with overdensity 
70 > A > 10 - 20 (the upper limit 70 is set by the local 
overdensity at the virial radius of a halo with a nonsin- 
gular isothermal density profile and core radius 0.1r V j r ). 
Using the distribution given by equation ((19(1 at z ~ 4—5, 
the clumping factor due to gas in regions with overden- 
sity between ~ 10 — 20 and 70 is 



(1 



^mii) 



70 



A 2 P(A)dA - 2 - 3. (23) 



HII 



15 - 16. 



(21) 



On the one hand, this value can only be taken as a lower 
limit since the numerical simulations used to obtain equa- 
tion ((19|) do not resolve halos with mass < 10 6 M Q and 
the dumpiness in these smaller scales may also contribute 
to Chii- On other hand, the clumping factor estimated 
above may be not significantly smaller than the reality, 
since most of the gas in halos with mass < 10 6 Mq, es- 
pecially in small halos, would be easily pho toevaporated 
ijBarkana fc Loeblll999l:IShapiro et alJl20Q4|) . and part of 
the gas in regions with overdensity between ~ 10 — 20 
and 70 may actually be ionized by photons from stars. 

Note that the contribution from equation ((23(1 does 
not include the contribution due to minihalos with mass 
between M rcs and M m j n , which are also parts of com- 
ponent iii. W ith a simple model of the photoevapo- 
ration process, lHaiman et alJ l|2001|) argued that these 
minihalos are potentially important sinks of ionizing 
photons. Using numerical simulations o f gas dynam- 
ics an d radiation transfer in minihalos, iShapiro et alJ 
(2004) argued that at z = 9 it may take only about 
(1 — 1.5) x 10 8 yr to photoevaporate these minihalos by 
a source at a distance of 1 Mpc away and emitting ion- 
izing photons at a rate of 10 56 s _1 (or equivalently, at a 
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distance of 10 kpc away and with emission rate 10 52 s" 1 ), 
and thus the photoevaporation process may reduce the 
required number of ionizing photons to reionize these 
minihalos. Our calculations in § 13.31 show that the av- 
erage ionizing photons emitted from stars or galaxies 
at a distance of 1 Mpc from those minihalos is about 
(0.5 - 1.1) x 10 57 (1 Mpc/r) 2 s^ 1 ~ (0.2 - 0.4) x 10 56 s" 1 
for model a and (0.8 - 1.6) x 10 57 (1 Mpc/r) 2 s" 1 ~ 
(0.3 — 0.6) x 10 56 s _1 for model b, respectively. These 
ra tes are si gnificant ly l ower than the required rate given 
bv lShaoiro et alJ l)2004|) to evaporate the gas in minihalos 
within a period of (1 — 1.5) x 10 8 yr. Therefore before the 
nuclear activity of the QSO turns on, despite that much 
of the gas in minihalos may have already been ionized 
and pushed out of the halo gravitational potential wells 
by the ionizing photons from stars, there should still re- 
main a considerable fraction of gas in minihalos. After 
the nuclear activity turns on, the remaining gas will be 
ionized by QSO photons. As will be discussed in § 0] 
the time during which the QSO has radiated around the 
observed luminosity iVp hs q before the cosmic time t(zo) 
is probably at most a few times 10 7 yr, which is shorter 
than the period (1 — 1.5) x 10 8 yr above; thus not all of 
the gas is evaporated at z ~ zq. Using the static approx- 
imation of the photoev aporation process in minihalos, 
IBarkana fc Loebl l)1999j) show that a significant fraction 
of the gas (about 20%-40% for a typical QSO spectrum 
at redshift z ~ 8 - 20) of big mini-halos (> 1O 6 M ) 
is still gravitationally bounded to the potential well af- 
ter the reionization of the universe is completed. Based 
on the tendency of the fraction with redshift shown in 
IBarkana fc Loebl Jj.999), this fraction could be larger at 
a lower redshift zq ~ 6.2 — 6.4. Below we simply assume 
that half of the gas still remains in minihalos and obtain 
the contribution of the minihalos to the clumping factor 
as follows: 
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(24) 

Our calculation shows that within a proper scale of tq = 
5 Mpc, the number of halos with mass between M res and 
Mmin passed by QSO photons at a random line of sight is 
on average 2 — 3 before the photons are all absorbed, so 
that the average clumping factor obtained above can be 
applied to the specific observed QSOs. Our calculation 
also shows that the fraction of hydrogen remaining in the 
minihalos is 
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dM ~ 3%. (25) 



Combining the contributions from components ii and 
iii, the clumping factor in the spheres considered in this 
paper (with tq = 5 Mpc at z ~ 6) is 

Chii = C HII + Chii ~ (34 — 38). 

To more accurately estimate the clumping factor would 
require very high resolution numerical simulations with 
sufficiently large volume and incorporating dynamical 
effects of photo-evaporation and radiation transfer and 
simultaneous consideration of the ionizing effects from 
both stars and QSOs. The clumping factor Chii may 



slowly decrease with time elapse because of the evolu- 
tion of the gas density distribution (e.g., due to photo- 
evaporation of some minihalos). Also, the distribution 
of (mini)halos within the Stromgren sphere is likely to 
be nonuniformly clustered and could further introduce 
fluctuations in Chii for different lines of sight. In § 0] 
we shall use not only the estimate above but also some 
other different values of the clumping factor to illustrate 
the results. 

4. THE EVOLUTION OF THE STROMGREN SPHERES 
AROUND THE HIGHEST REDSHIFT QSOS 

4.1. Models 

After the nuclear activity of the highest redshift QSO 
turns on, the neutral hydrogen remaining in high-density 
regions not ionized by stars will be further ionized by 
QSO photons and thus the neutral hydrogen fraction in 
the IGM surrounding the QSO will decrease. A spherical 
ionization front centering on the QSO will appear, sepa- 
rating the inside highly ionized HII region (due to QSO 
photons) and the outside partly ionized region (due to 
stars). The expansion of the apparent radius of the ion- 
ization front detected by observers (or the Stromgren ra- 
dius; see the apparent shape of the Stromgren sphere in 
lYull2004|) . thii, can be described by the following equa- 
tion for the number of hydrogen atoms ionized by the 
QSO in the Stromgren sphere: 



4-7T d(x H i(n H )r HII ) 



dr 

4-7T 



aBC H n(riH) r mi , 



(26) 



where t = t — ti is the time that the QSO has passed at 
cosmic time t since its nuclear activity was triggered at 
cosmic time U, N v h s ,Q{T) is the evolution of the QSO ion- 
izing photon emission rate, and «b = 2.6 x 10 -13 cm s 
is the hydrogen recombination coefficient to excited lev- 
els of hydroge n at temperature T = 10 4 K (see Donahue 
fc ShuU 1987: iShapiro fc Giromd Il987t ICen fc Haimail 
120001 iMadau fc Reesll2000HYull2004|) . The terms on the 
right-hand side of equation l)26|) account for the ioniza- 
tion due to QSO ionizing photons and the recombination 
within the Stromgren sphere, respectively. The role of 
the recombination can be characterized by a recombina- 
tion timescale defined as follows: 



Tree = x m (Chii (tie.) a B ) 1 

Xm / 30 



4.8 x 10 7 yr 



(i + <W) VCi 



'HII 



7.42 

1 + z 



(27) 



The recombination is negligible if r <C T rcc , while impor- 
tant if r > Tree- 
In equation 126|) the parameters xhi, (^h)j and Chii 
are the average values within the dense region around a 
QSO with radius rnn at cosmic time t, and their changes 
during the time interval rnn/c of the QSO photon prop- 
agation in the Stromgren sphere are assumed to be neg- 
ligible. Their changes caused by the Hubble expansion 
are also negligible. For simplicity, we shall also neglect 
their dependence on rnn in the calculations below and 
set them to be the values averaged within a fixed scale, 
e.g., ran = rq at zq, as analyzed in §[3] This simplifica- 
tion will not affect our conclusions. 
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For the evolution of the QSO ionizing photon emis- 
sion rate A p h s ,Q, below we assume two simple mod- 
els: in model i A p h s ,Q = Ap hs Q i s constant, and in 
model ii A p h s ,Q increases with time exponentially, i.e., 
iVphs.qM = iV p hs , Q exp[(T - t q )/t s ] (0 < r < tq), 
where tq = t(zq) — U is the age of the QSO at zq, 
ts — 4.5 x 10 7 yr[e/0.1(l — e)] is the Salpeter timescale, 
and e is the m ass-to-energy convers i on efficiency and 
is ~ 0.1 (e.g.. IGammie et al.l 12001 lYu fc Lul I2004at 
lYu fc Tremainel2002fl . The ionizing photon emission rate 
is proportional to the QSO luminosity if thei r intrinsic 
continu um spectra do not differ significantly. lYu fc Lul 
( 2004a) show that the exponential increase of the QSO 
luminosity with a characteristic Salpeter timescale (that 
is, QSOs radiate at a rate close to their Eddington lu- 
minosity) appears to be generally consistent with the ex- 
pected demographic relations between the QSO luminos- 
ity function and massive dormant black holes in nearby 
galaxies. It is likely that the luminosity evolution of a 
QSO also experiences a declining-luminosity phase (e.g, 
with the consumption of the accreting materials at the 
end of the nuclear activity), but this phase generally does 
not dominate the main population of luminous QSOs and 
is not considered in this paper, for simplicity l)Yu fc Lul 
l2004a|h Models i and ii are almost the same if tq <ts. 

The highest redshift QSOs are also observationally sug- 
gested to be accreting at a rate Mbh close to the Edding- 
ton limit ( Wil lott et al.1 120031) . With this observational 
constraint, in model i, the QSO age tq should be not 

longer than Mbh/ Mbh ~ ts (for this reason we extend 
the TQ-axis in Fig. at most to a few times 10 7 yr, 
rather than to 10 s yr shown in Fig."^), and the cumula- 
tive ionizing photons emitted from the central QSO are 
at most 1-2 per hydrogen atom in the IGM within their 
Stromgren spheres (if e ~ 0.1), which is significantly less 
than those from stars (see S I3.3[I . Although in model ii tq 
may be much longer than ts, before t(zQ) most ionizing 
photons are emitted within the period from £(zq) — t$ 
to t{zQ). The time interval ts corresponds to a redshift 
interval of Sz ~ 0.3 at zq ~ 6.2 — 6.4. Cumulatively, 
most ionizing photons from stars before t(zo) are emit- 
ted before t(zQ + Sz), and before t(zQ +5z) most ionizing 
photons are not from QSOs but stars (see Fig. [""J, which 
helps to justify the isolation of the ionization effect due 
to stars and that due to QSOs done in this paper (i.e., 
first consider stars and then QSOs). 

Depending on whether the QSO shining time tq is 
much shorter than the recombination timescale r rec (and 

the Salpeter timescale ts for model ii of A p h s ,Q), equa- 
tion l|26|) and its solution can be simplified into the fol- 
lowing two cases: 

• Case (1): if tq <C r rec (and tq <§; ts for model 

ii of iV p hs,Q), the recombination term in equation 
(|26|) is negligible and the apparent size of the highly 
ionized HII region at r = tq is given by 



' HII 



3/ TQ ^ P hs,Q(r)dr 



AzT<X\ 



^^ p hs , Q TQ ' 
47ra;Hi(n-H) 



(28) 



(29) 



In this case, given Np hs Q and (uh), the size of 

r Hii ( t q/ 2; hi) 1 ^ 3 depends on the QSO age tq as 
well as the neutral hydrogen fraction xhi- There is 
little difference in the results from models i and ii. 



Case (2): if tq 3> T rcc , the recombination within 
the highly ionized HII region is approximately bal- 
anced by the emission of the ionizing photons from 
the QSO, and we have 



'3JV° 



'HII- 



phs,Q T re 



1/3 



47TXHi(n-H) 



3/V° 



1/3 



I 47ra B CHii(nH) 



rs 



for model i and 



r H n-(— ^-J - 

T rec + TS 



1/3 



(30) 



(31) 



for model ii. The definition of rs above is similar 
to the classical formula of the Stromgren radius in 
stellar astronomy. In this case, given Np^Q and 
does not depend on the exact value of 



\ n H) 



HII 



the QSO age tq. If ts S> T rec in model ii, there is 
little difference in the results of the two models and 
r Hii( cx Cffl/ 3 ) increases with decreasing Chii- The 
dependence of Chii on xhi is probably weak since 
most hydrogen atoms are actually located in low- 
density regions, while Chii is mainly determined 
by the hydrogen i n high-density regions (as illus- 
trated in Fig. 2 in iMiralda-Escude et ai]l2000fl . If 
ts <C T rcc in model ii, obtained from model ii is 

l/3\ 



h 'HII 

smaller than that from model i, and r^ u ((x x^ 6 ) 



depends mainly on xhi in model ii instead of Chii 
in model i. 

Below using the A p h s ,Q given by models i and ii, we 
solve equation (|26|l to illustrate the general properties 
of the evolution of rnn by artificially setting different 
values for the clumping factor Chii and the neutral hy- 
drogen fraction xhi ■ The solutions of r HII are shown as 
a function of the assumed QSO age tq in Figure[21 Pan- 
els (a) and (b) give the results for models i and ii of 
-Wphs,Qj respectively. In both panels, the QSOs are as- 
sumed to be at zq = 6.3 and have N® hs q = 10 57 erg s _1 , 
and we set <5 rml ~ 0.25, Chii=10, 30, and 50 (dot- 
ted, dashed, and solid lines, respectively), and x^i = 
10" 2 , 10" 15 , 10" 1 , 10~ - 5 , and 1 (top to bottom for each 
line type). As seen from both panels, at the short-age 
end (e.g., < 10 7 yr), r^jj increases with increasing tq 
and decreasing xhi as described by case 1 above in which 
the effect of recombination is insignificant (see eq. I29fl : 
while at the long-age end, with increasing QSO ages, r HU 
roughly approaches a constant and the constant increases 
with decreasing Chii as described by case 2. Given Chii, 
the value of the constant is generally independent of xhi 
but may decrease with increasing xhi if £hi is sufficiently 
large (close to 1 so that ts < T roc ) in Figure 2b for model 
ii (see eas. 13*01 and l3*T|) . 
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Fig. 2. — The expected apparent size of the Stromgren sphere 
r^jj as a function of the QSO age tq. The QSOs are assumed to be 
at redshift zq = 6.3. Different line types represent different clump- 
ing factors used (dotted, dashed, and solid lines for Chii =10, 30, 
and 50, respectively). For each line type, we set the neutral hydro- 
gen fraction x m = 10 -2 , 10" 1 ' 5 , 10 _1 , 10 -0 5 , and 1 from top to 
bottom. Panels (a) and (b) are for different models of the evolution 
of the QSO ionizing photon emission rate N p ^ s Q (models i and ii 
in § 14.11 respectively). As seen from this figure, if tq is signifi- 
cantly short, r^jj increases with increasing tq and decreasing xjji 
(see eg. 1271 . If tq is long enough (e.g., > 10 7 yr), r^ n approaches 
a constant, independent of the detailed values of tq. The constant 
is insensitive to the exact values of xjji and the evolution models 
of ^Vpha.Q used and is mainly determined by the clumping factor 
(the highest redshift QSOs considered in this paper belongs to t his 
case ), u nless xjji is big e noug h (close to 1) in panel (b) (see eqs. 1301 
and||n}. See details in S lO 



4.2. Comparison with observations 

For the highest redshift QSOs considered in this paper 
(<5 rmi ~ S rQ ~ 0.25, < 0.1-0.2, Chii ~ 36, z Q -6.2- 
6.4; see details in § 13.31 and § I3.4f> . the averaged recom- 
bination timescale of the hydrogen in their surrounding 
regions is about r rcc <5x 10 6 yr. Note that the average 
lifetime of the main population of QSOs (with comoving 
number density peaked at z —2-3) has been determined 
to be mfe > 4 x 10 7 yr <|Yu fe Lull2004at lYu fe Tremalnd" 
120021 see also Martini 2004 for a review of the QSO life- 
time). If the lifetimes of the highest redshift QSOs are 
also > 4 x 10 7 yr, the probability that their ages are 
shorter than r rcc is therefore only — T rcc /riif < 10%, 
and the ages of the majority of the observed highest red- 
shift QSOs are more likely to be around or longer than a 
few times 10 7 yr, which is significantly longer than r rec . 

Even without the constraints of the QSO lifetime from 
other methods and an accurate estimate of Chii, by es- 
tablishing a statistical method relating the distribution 
of the observed Stromgren radii with the probability dis- 
tribution of the ages of the obs erved QSOs (som ewhat 
similar to the manipulations in lYu fe Lul I2004albf) . one 
can still check or rule out the possibility of whether 
the lifetime of the highest redshift QSOs is shorter 
than T rcc or whether the solution of r HII belongs to 
case 1. For example, the observational distribution of 

((nH)/^Vphs,Q)~ 1 ^ 3r Hii.obs can be estimated from a sam- 
ple of QSOs having Gunn-Peterson troughs. This distri- 
bution should be consistent with a constant if most of the 
r Hii solutions belong to case 2, and consistent with the 

1/3 

distribution of Tq obtained by a random choice of tq 



for each QSO if most of the r HII solutions belong to case 
1 and xm is irrelevant to tq and does not differ much for 
different QSOs in the sample. Our preliminary results of 
applying this statistical method to the three highest red- 
shift QSOs with measured Stromgren radii (in Tab. [T]) 
have shown that the possibility of Tuf c < T roc is less than 
15%. For simplicity, the details of this method are not 
presented in this paper and will be deferred to future 
work with an increasing number of the detected highest 
redshift QSOs having Gunn-Peterson troughs. 

We conclude here that for most highest redshift QSOs 
the solution of r HII belongs to case 2 and is insensitive to 
the exact values of the QSO age. For the highest redshift 
QSOs, we also have r rcc <C t$ (in model ii), and thus r HII 
is also insensitive to the neutral hydrogen fraction and 
the detailed evolution models of -/V p hs,Q (see eas.lBTjland 
I31|) . The r HII value is mainly determined by the QSO 
ionizing photon emission rate N® hs q and the clumping 
factor Chii- Below we simply show the result obtained 
from model i. 

We plot r HII as a function of the ionizing photon emis- 
sion rate of QSOs -/Vp hs in FigureEt- The observation- 
ally determined Stromgren radii of three QSOs and their 
ionizing photon emission rates (see Tab. 1) are shown 
by filled circles. The model results obtained by setting 
S rHU ~ S rQ ~ 0.25, Chii = 36, and sehi ^ 0.1 — 0.2 
are shown by the dotted and solid curves, and by an 
open circle (in the range of r HII — 4.9 — 5.9 Mpc) for the 
source SDSS J1048+4637, which does not have a mea- 
sured Stromgren radius in the literature. The different 
solid lines are for QSOs at different redshifts (zq = 6.2, 
6.3, and 6.4; top to bottom). We also illustrate the effects 
of Chii by showing the results obtained by three other 
values of Chii=20, 30, and 45 in Figure OK (top, middle, 
and bottom dotted lines, respectively; zq = 6.3). The 
relation between the model and observational results for 
the three QSOs are more explicitly shown in Figure \S]p 
(see solid circles). As seen from Figure the model re- 
sults are remarkably consistent with observations, which 
supports that the lifetimes of the highest redshift QSOs 
are around or longer than a few times 10 7 yr, as is the 
lifetim e of the main populati on of QSOs (e.g.. lYu fc LiJ 
I2004al lYu fc Tremaind 12002(1 . The remarkable consis- 
tency also supports that the observational sizes of the 
Stromgren spheres of the highest redshift QSOs can be 
self-consistently explained by the fact that they are lo- 
cated in rare overdense regions and that the neutral hy- 
drogen fraction in these regions is in the range from sev- 
eral percent to 10%-20% (before considering the reion- 
ization due to QSO photons). A much smaller neutral 
hydrogen fraction (e.g., <C 3%) will not be consistent 
with the clumping factor used, since Chii ~ 36 implies 
that a significant fraction of the neutral hydrogen (see 
cq . 125(1 remains in (mini)halos (which was neither ionized 
by the ionizing photons from stars nor photo-evaporated 
out of the halos) and contributes to Chii (see ea. l25ll . Un- 
less Chii is much smaller than the estimate in this paper 
(e.g., because of stronger photo-evaporation of baryons 
in minihalos), a much shorter QSO lifetime (<C 10 7 yr) 
will not be consistent with the observations of r ml obs 

since r HIT decreases with decreasing tq at short-TQ ends 
(see Fig. EJ). 
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Fig. 3. — Panel (a): The size of the highly ionized HII regions 
r Hll °f * ne highest redshift QSOs as a function of the ionizing 
photon emission rate Np hs q. The filled circles represent the ob- 
servational data of the three QSOs with measured r^ n obs given in 
the literature (see Tab.0. The open circle gives our model results 
of r^jj mod for the QSO J1048+4637 obtained by solving equation 
1261 and setting N® hs q to be the mean of its upper and lower 
limits (listed in Table 0, and its estimated error bar in r^ n mod 
represents the range of the upper and lower limits of N® ha q (simi- 
larly for the model results of other sources shown in panel b) . The 
dotted and solid curves also represent the model results. The dot- 
ted lines give the results for QSOs at zq = 6.3, and the solid lines 
give the results for QSOs at zq = 6.2, 6.3, 6.4 from top to bot- 
tom, respectively. The top and bottom dotted lines are obtained 
by setting Chii=20, 30 and 45, respectively; and the solid lines 
and the open circle are obtained by setting Chii = 36. Note that 
for low-luminosity QSOs (with N® hs q smaller than a few times 
10 56 s — 1 ) located in less dense regions, Chii = 36 could be an 
overestimate and the range of observed ?* HII could be higher than 
those shown by the dotted and dashed lines and up to the dashed 
line (obtained by Chii = 20; see the region labeled by "?"). Panel 
(b): the model results of the apparent size of the Stromgren sphere 
r HH mod versus the size obtained from observations for the highest 
redshift QSOs r Hn obs (solid circles). The solid line is a refer- 
ence line. The remarkable consistency of the model results with 
observations supports that the observed Stromgren radii can be 
self-consistently explained by the dense environment around the 
highest redshift QSOs, the corresponding biased structure forma- 
tion, star formation, and reionization processes analyzed in this 
paper, and the QSO lifetime longer than a few times 10 7 yr. If the 
observational r^ TI obs is systematically underestimated, then the 
realistic Chii ma y De lower than ~ 36 and the photo-evaporation 
process in minihalos may be stronger than that assumed in this 
paper. See details in S fH 



The solution of r HII in model ii is lower than that in 
model i by at most 1 - [r s /(r rcc + t s )} 1/3 ~ t icc /(3t s ). 
An independent constraint on the timescale t$ may be 
obtained by comparing the result of model ii with the 
observation, as shown in Figure |3Jd, but this may require 
precise measurements of the observational Stromgren ra- 
dius r^ II obs . The current observational error of ^ II obs 

(^ r Hii obs ~5%-30%; see Tab. and the required con- 
sistency with the observation may give the constraint to 
be r s /r roc > l/(35r° IIi0bB ) -1-7. 

4.3. Discussions 
The model results of r^ u mod above are obtained by 
using the average clumping factor (over all the lines of 
sight). A significant contribution to the clumping fac- 
tor comes from the neutral hydrogen in halos. In sparse 
regions, the ionizing photons may have been absorbed 
before they encounter any halo with mass greater than 
M m i n at a significant fraction of the lines of sight. Thus 
the average contribution of component iii to Chii may 
be an overestimate when applying it to a specific line 
of sight (see discussions in § 13. 4|) . Hence, the model re- 
sults above are usually appropriate to be applied to the 
observed high-luminosity QSOs located in dense regions, 
but are underestimates for low-luminosity QSOs (e.g., 
ionizing photon emission rate ~ 10 56 s _1 at zq —6.2- 
6.4) located in less dense regions. The expected range 
of 

r Hii f° r low-luminosity QSOs can be estimated in the 
following way. By starting with their number density 
Nq (obtained by extrapolating the luminosity function 
of the highest re dshift QSOs with a po wer-law index of 
3.2; see details in lFan et al.ll2003ll2004j) at the faint end 
and doing similar quantitative analysis for the environ- 
ment around the low-luminosity QSOs as done for the 
high-luminosity QSOs above, we find that the clumping 
factor Chii is around 40 for those lines of sight along 
which one halo with mass greater than M m i n can be 
encountered within the Stromgren sphere, and around 
20 for other lines of sight (e.g., 50% of all the lines of 
sight) along which no halo with mass greater than M m i n 
is encountered (the details of the calculations are not 
given here for simplicity). Therefore, at a random line of 
sight, the size of the Stromgren spheres around the low- 
luminosity QSOs are roughly in the range of the results 
obtained by using Chii ~ 20 (dashed line in Fig. [3K) an d 
Chii — 36 (see the region labeled by "?" between the 
dashed line and the solid lines in Fig. |3Ji) . We believe 
that comparison of the range with future observations of 
the Stromgren radii of low-luminosity QSOs would fur- 
ther give some tests or constraints on the analysis or as- 
sumptions made in this paper, for example, on whether 
the major part of the clumping factor is due to the re- 
maining neutral hydrogen in halos. 

In addition, we note that a higher neutral hydrogen 
fraction (xht > 0.3 if tq > 10 7 yr) is obtained in 
iWvithe k Loebl l)2004a|) . but it is obtained by using a 
much smaller clumping factor and without considering 
the detailed effects of the structure and star formation in 
the dense environment surrounding the highest redshift 
QSOs. Also, as mentioned in § |H th e Stromgren radius 
of SP SS 1030+0524 estimated by IMesineer fe HaimarJ 
(2004) through the LyS trough in the QSO spectra 
(— 6 Mpc) is larger than the value (~ 4.6 Mpc) listed 
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in Table ^ A possibly larger observational Stromgren 
radius may suggest that the realistic clumping factor for 
this QSO is smaller than the estimate in this paper (e.g., 
Chii may be as low as ~ 20; see Fig. |2J. A smaller Chii 
may be caused by a stronger photoevaporation process 
in minihalos than that assumed in this paper. If more 
minihalos were photoevaporated, besides Chii, the neu- 
tral hydrogen fraction to be ionized by the central QSO, 
xhi, may also decrease. In short, most of the analysis 
in this paper will not be qualitatively affected even if 
the Stromgren radii estimated from current observations 
deviate slightly from reality. 

As indicated in Figure ^ m the rare overdense regions 
around the highest redshift QSOs, the ionizing photon 
production per hydrogen atom in the IGM at z ~ 7 is 
already comparable to or higher than that of the cosmic 
average at z ~ 6 (when the reionization is complete). 
However, the reionization in the rare overdense regions 
around the highest redshift QSOs is not complete at a 
redshift significantly higher than zq ~ 6.2—6.4 since neu- 
tral hydrogen must remain in the high-density regions to 
be ionized by QSO photons to satisfy the observational 
constraint from the Stromgren radius analyzed in this pa- 
per. The incomplete reionization in the overdense regions 
is partly caused by the fact that they have more high- 
density gas and that the gas is easier to recombine. This 
environmental effect has also been revea led in the nu- 
merical simulation bv iCiardi et alJ l)2003f) . which shows 
that the reionization process in overdense regions may 
be complete at a time later than that in cosmic average 
regions, although the ionizing photon production per hy- 
drogen atom in the IGM is higher in overdense environ- 
ments than in cosmic average re gions. According to this 
environmental effect revealed bv lCiardi et al.1 l)20f)3fl . the 
cosmic average neutral hydrogen fraction at z ~ 6.2 — 6.4 
should be not higher than 10%-20% (the upper limit of 
Xhi surrounding the highest redshift QSOs). The obser- 
vational flux upper limit of the Gunn-Peterson troughs 
has constrained the lower limit of the neutral hydrogen 
fraction (in mass ) of the universe at z ~ 6.2— 6. 4 to be 1% 
l|Fan et al.H2002lh Therefore, the cosmic average neutral 
hydrogen fraction at z ~ 6.2 — 6.4 may be only a few per- 
cent. This low neutral hydrogen fraction alleviates the 
apparent conflict in reionization betwee n the constraint 
from t he highest redshift QSOs (e.g., iWvithe fc Loebl 
l2004a|) and that from the polarization spectrum of the 
cosmic microwa ve background (i.e. . an early reionizatio n 
at z ~ 15 or so: iKogut et"aDl2r)03t ISnergel et al.ll2003ft . 

5. CONCLUSIONS 

In this paper we have investigated the dense environ- 
ment, the reionization process, and the evolution of the 
Stromgren sphere around the highest redshift QSOs hav- 



Abel, T., Anninos, P., & Norman, M. L. 1997, New A, 2, 181 
Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 
Bajtlik, S., Duncan, R. O, & Ostriker, J. P. 1988, ApJ, 327, 570 
Barkana, R. 2004, MNRAS, 347, 59 
Barkana, R., & Loeb, A. 1999, ApJ, 523, 54 
Barkana, R., & Loeb, A. 2001, Phys. Rep., 349, 125 
Barkana, R., & Loeb, A. 2004, ApJ, 609, 474 
Becker, R. H., et al. 2001, AJ, 122, 2850 

Benson, A. J., Nusser, A., Sugiyama, N., & Lacey, C. G. 2001, 
MNRAS, 320, 153 



ing Gunn-Peterson troughs [z > 6.1), and we have pro- 
vided constraints on the intrinsic properties of QSOs and 
the reionization history of the universe by comparing the 
observed Stromgren radii with model results. 

We have shown that the structure formation and con- 
sequently the gas distribution and star formation in the 
overdense regions around the highest redshift QSOs are 
biased from the cosmic average. Before the nuclear activ- 
ity of the QSO turns on and its ionizing photon emission 
rate is high enough, the reionization in the overdense re- 
gion surrounding the QSO is mainly contributed by stars, 
and starts from relatively low density subregions. Using 
some simple models of star formation and some analysis 
on the photoionization process inside and outside of halos 
(including minihalos), we have argued that a significant 
fraction of hydrogen in the Stromgren sphere around the 
QSOs is ionized by the ionizing photons from stars, and 
only about several percent to at most 10%-20% of hydro- 
gen is left (e.g., in minihalos, halos, or high-density sub- 
regions) to be ionized by the QSO photons. The cosmic 
average neutral hydrogen fraction at z ~ 6.2 — 6.4 is also 
be smaller than the upper limit of 10%-20% obtained 
for the overdense regions and may be only a few per- 
cent, since overdense regions contai n more high-densit y 
gas and are more difficult to reionize (|Ciardi et alJ2003|) . 
We have analyzed the clumping property of hydrogen 
ionized by QSOs and studied the evolution of the appar- 
ent size of the Stromgren sphere. We have found that 
if the QSO lifetime is about or longer than a few times 
10 7 yr, as is the lifetime of the main population of QSOs 
(with comoving number density peaked at z ~2-3), the 
expected Stromgren radii from our models are very con- 
sistent with observations. With such a QSO lifetime, 
the ages of most of observed QSOs are long enough that 
the QSO photon emission is balanced by the recombina- 
tion of the hydrogen ionized by QSOs in their Stromgren 
spheres, and the expected Stromgren radii from the bal- 
ance are independent of the detailed values of the QSO 
ages. We also point out a statistical method involving 
a larger sample of highest redshift QSOs having Gunn- 
Peterson troughs in future observations, which may po- 
tentially check or rule out the possibility that the QSOs 
have a shorter lifetime (e.g., < 10 7 yr) even without an 
accurate estimate of the hydrogen clumping property. 
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Table 1. Sample of QSOs that have the Gunn-Peterson trough 



Name 


Z Q 


Afl450 


JVP, n 

phs,Q 

(lO^s- 1 ) 


ZHII.obs 
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„0 

' Hll.mod 










(Mpc) 


(Mpc) 
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Table 2. Mass fraction of baryons in formed stars and ionizing photon 
emission rates from stars 



Model 


(xicr 2 ) 


9star(2Q) 

(xlO" 2 ) 


2 

9star(z) = 9%S( Z Q) 


N h l aB . (zq) 

pns,star V 

(10" s" 1 ) 


N phs,star(2 ~ 5) 

(lO^s" 1 ) 


(a) 
(b) 


1.3-1.4 
1.0-1.2 


0.8 
0.6 


4.5-4.8 
5.0-5.3 


0.5-1.1 
0.8-1.6 


0.3-0.6 
0.4-0.8 



Note. — Models (a) and (b) are two different models for the star formation efficiency (i.e., the mass fraction of baryons to form stars) 
in halos (see details in § 13.31 . The g s t ar (zq) gives the mass fraction of baryons that are in formed stars in the rare overdense regions with 
proper radius tq = 5Mpc around the highest redshift QSOs, where the QSO redshift 2q ~ 6.2 — 6.4; and the g s tar(2Q_) gives the cosmic 
mean fraction at the same redshift. The fourth column gives the redshift at which g s tar(z) = Sstar ( 2 q) (indicated by the arrows in FigureHl 
The redshift is close to z ~ 5 for both models. The ^fp^ B staT ( z Q) gives the ionizing photon emission rates from stars in the overdense 
regions associated with the highest redshift QSOs at redshift zq; and the Af phs star (z ~ 5) gives the cosmic mean rate at z ~ 5. We have 
^ph" tar(*Q) - 2A r phs , s tar(2 ~ 5) for both models. See also Fig. [T] and §EJ 



